Heavy quark energy loss and D-Mesons in RHIC and LHC energies 
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We obtain the radiative energy loss of a heavy quark in a deconfined medium due to radiation 
of gluons off them using a recently derived generalised gluon emission spectrum. We find that the 
heavy flavour loses energy almost in a similar fashion like light quarks through this process. With 
this, we further analyse the nuclear modification factor for D-meson at LHC and RHIC energies. 
In particular, the obtained result is found to be in close agreement with the most recent data from 
ALICE collaboration at 2.76 ATeV Pb-Pb collisions. We also discuss the nuclear modification factor 
due to the collisional energy loss. Furthermore, the result of non-photonic single electron from the 
decay of both D and B mesons is compared with the RHIC data at 200 AGeV Au-Au collisions, 
which is also in close agreement. 



INTRODUCTION 

The purpose of ongoing relativistic heavy ion collisions 
is to understand the properties of nuclear or hadronic 
matter at extreme conditions. A primary aim lies in the 
detection of a new state of matter formed in these colli- 
sions, the quark-gluon plasma (QGP), where the quarks 
and gluons are liberated from the nucleons and move 
freely over an extended region rather than over a lim- 
ited hadronic volume. Various diagnostic measurements 
taken at CERN Super Proton Synchrotron (SPS) [lj in 
the past and at BNL Relativistic Heavy Ion Collider 
(RHIC) [2tS| in the recent past have provided strong 
hints for the formation of QGP within a first few fm/c 
of the collisions through the manifestation of hadronic 
final states. New data from heavy-ion experiments at 
CERN Large Hadron Collider (LHC) 0-[ll| bave further 
indicated the formation of such a state of matter. 

One of the important features of the plasma produced 
in heavy-ion collisions is suppressed production of high 
energy hadrons compared to the case of pp collisions, 
called jet quenching. The term 'jet quenching', gener- 
ally, ascribes to the modification of an energetic parton 
due to its interaction with the coloured medium while 
passing through it. The basic idea is that the scales of 
hard (high-pjj processes and the medium interactions in 
the context of heavy-ion collisions, are very distinct in 
accordance with the uncertainty principle. This provides 
the fact that the high-pj_ parton production in A— A colli- 
sions can be computed using perturbative QCD (pQCD), 
which is quite close to the vacuum rate scaled for binary 
N — N collisions in an A — A collision. The effect of 
medium is then treated as a final state interaction which 
is taken into account through the modification of the 
outgoing parton fragmentation pattern due to parton- 
mcdium interactions. 

The heavy-ion program at BNL RHIC |8j has clearly 
revealed that the phenomenon of jet quenching is mainly 
caused due to the energy loss of the initial hard parton via 



collisional and radiative processes, prior to hadronisation. 
The indication for jet quenching in heavy-ion program 
at CERN LHC has also been |9Hll| observed recently. 
The energy loss encountered by an energetic-parton in 
a QCD medium reveals the dynamical properties of that 
medium and presently is a field of high interest in view of 
jet quenching of high energy partons; both light 12H28J 



and heavy quarks [16, 29-51]. Naively, one imagines that 



the amount of quenching for heavy flavours jet should be 
smaller than that of light flavours due to the large mass 
of heavy quarks. However, the single electron data at 
RHIC [8] exhibit almost a similar suppression for heavy 
flavored hadrons compared to that for light hadrons. 

A first attempt to estimate the radiative energy loss of 
heavy flavours in a QGP medium was made in Ref. 32 1 
by using the Gunion Bertsch (GB) formula [52j of gluon 
emission for light quark scattering and appropriately 
modifying the relevant kinematics for heavy quarks. 
Later the GB-like formula for heavy quarks was recon- 
sidered in Ref. 33j by introducing the mass in the matrix 
element |34| but only within the small angle approxima- 
tion. Due to this mass effect, a suppression, known as 
'dead cone' effect, in the soft gluon emission off a heavy 
quark was predicted in comparison to that from a light 
quark. This resulted in a reduction of heavy quark energy 
loss induced by the medium [33j , which is limited only to 
the forward direction. However, such a gluon radiation 
spectrum with a dead cone factor, only applicable to the 
forward direction, was also used in the literature [35j, |37| 
uniformly for the full range of the emission angle (i.e., 
both forward and backward direction) of gluon to cal- 
culate the heavy quark energy loss in the medium. This 
can lead to a unphysical result at large angle radiation, as 
discussed as well as shown in Re f. I53 j|. Further attempts 
were also made in the literature HJH, H, EH 0, 0, U5l| 
to improve the calculation of heavy quark energy loss 
with various ingredients as well as restrictions. In some 
cases the energy loss for charm quark was found to be dif- 
ferent than the light quark. The subject of heavy quark 



energy loss is not yet a settled issue and requires more 
detailed analysis. 

In a very recent work 53] the probability of gluon emis- 
sion off a heavy quark has been generalised by relaxing 
some of the constraints, e.g., the gluon emission angle and 
the scaled mass of the heavy quark with its energy, which 
were imposed in earlier calculations [33J, |34| . It resulted 
in a very compact and elegant expression for the gluon 
radiation spectrum off a heavy quark (e.g., Qq — > Qqg) 
as 53] 
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where the transverse momentum of the emitted massless 
gluon is related to its energy by k± = uj sin 0, and the 
rapidity, rj = — ln[tan(0/2)], is related to the emission 
angle, and the generalised dead cone is given by 
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Now, 



the Mandelstam variable s is given as, s = 2E + 
2E^E 2 - M 2 - M 2 , with E and M, respectively, the 
energy and mass of the heavy quark. C 'a is the Casimir 
factor for adjoint representation and a s is the strong 
coupling constant. In the small angle limit, 9 <C #o(= 
M/E) <C 1, the dead cone in (2) reduces to that in 
Ref. [33j, |34| as (1 + 2 /9 2 )~ 2 whereas for massless case 
it becomes unity and ([T]) reduces to the GB formula [52j. 
The gluon spectrum for the process, Qg — > Qgg, can also 
be found in Ref. [53|. We also note that the gluon emis- 
sion spectrum in (JXJ) is obtained in Feynman gauge. The 
same result is also obtained using light-cone gauge. 

In Fig. [TJ a Monte Carlo simulation of the above sup- 
pression factor (2) (i.e., the scaled gluon emission spec- 
trum off a heavy quark with that of light quark) is dis- 
played. It reveals a forward-backward asymmetry which 
encompasses the fact that the gluon emission off a heavy 
quark is as strong as that of light quark at the large an- 
gles (backward direction) whereas it is suppressed due to 
nonzero quark mass at the small angles (forward direc- 
tion). However, if the energy of the heavy quark is large 
compared to its mass, the effect of dead cone diminishes, 
both heavy and light quark are expected to lose energy 
almost similarly. This result can have important conse- 
quences for a better understanding of heavy flavour en- 
ergy loss in the context of heavy-ion collisions at RHIC 
and LHC. In this article we intend to use the gluon ra- 
diation spectrum in Ref. [53] to obtain the heavy flavour 
energy loss and attempt to understand the suppression 
of heavy flavoured hadrons in heavy-ion collisions. 
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FIG. 1. (color online) A Monte Carlo simulation for the sup- 
pression factor in © in the full domain of gluon emission 
angle, 9, off a heavy quark for the scaled mass m = M* — 0.3. 
This actually represents a two dimensional view of the scaled 
gluon emission probability off a heavy quark with that of a 
light quark as given in fl]). We consider the direction of prop- 
agation of a heavy quark is from left to right along the hori- 
zontal axis and collide with medium partons at the origin of 
a circle of unit radius. This simulation has been performed 
by throwing points at random directions within the full do- 
main of 6 but with a probabilistic weight T>(0), which would 
then correspond to a point randomly on the selected 6-line as 
a 'red plus' inside the circle of unit radius. The shade with 
red pluses represents the soft gluon emission zone whereas the 
conical white zone in the forward direction indicates a dead 
cone for gluon emission due to the mass effect. 



tive) scattering is undoubtedly the most important and 
interesting one. A number of different energy loss models 
has also been formulated in the literature (for review see 
Refs [27|, |28[). The basic differences among the different 
models are the various constraints (e.g., kinematic cuts, 
large angle radiation etc.) implemented to make the cal- 
culations manageable. In this section we define the rate 
of radiative energy loss of a parton with energy E, due to 
inelastic scatterings with the medium partons in a very 
canonical way as 



dE 

da: 



M 

(A)' 
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where (uj) and (A) are the mean energy of emitted gluons 
and the mean free path of the traversing quark, respec- 
tively. 

Among the set of variables [k±,f],uj] in ([1} any two 
together are sufficient to completely describe an emitted 
gluon. For convenience we now change the variable duo 
from [fcj.,77] to [u3,rj\ as 



RADIATIVE ENERGY LOSS 

Among the interactions that a charged particle under- 
goes, as it traverses a dense matter, inelastic (i.e. radia- 
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gluons from the spectrum as 
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Other important quantity in §3§ is the mean free path 
(A), which is the average distance covered by the travers- 
ing quark between two successive collision, followed by a 
soft gluon radiation. The magnitude of mean free path 
depends on the characteristics of the system in which the 
energetic particle is traversing, and it is defined as 



(A) = l/(cr 2 -> 3 Pqgp) 



(6) 



where a 2 -^3Pq eP = Pq< 7 Q q (q)^Q q (q)g + PgO'Qg^Q gg , o"2->3 is 
the cross section of relevant 2 — ► 3 processes and p qgp is 
the density of QGP medium which acts as a background 
containing target partons, for the high energetic projec- 
tile quark. We also note that the Landau-Pomeranchuk- 
Migdal (LPM) interference correction may be marginal, 
which we would estimate below based on the formation 
time of the emitted gluon along with the kinematical re- 
strictions. Now, we recall the total cross section for 2 — >• 3 
processes as given in Ref. [54J as 
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= 4 C A a\ J -^ dqi J i duj Jv d V , (7) 

where q±_ is the transverse momentum of the exchanged 
gluon. Combining ([6]) and ([7|) the energy loss in ([3]) can 
be written as 
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where a factor of 2 has been introduced in 77 integral to 
cover both upper and lower hemisphere. We note that 
for T> — 1, © becomes equivalent to the massless case. 

At this point it is important to note that the hierarchy 
employed in obtaining ([T]) in Ref. [53| reads as 



(9) 



y/s,E^> V |*| ~9i»w>fci»mi) , 



where s, u, t are the usual Mandelstam variables and mr) 
is the Debye screening mass of the thermal gluons. Based 
on the above hierarchy we obtain the kinematic cuts ex- 
plicitly on energy-momentum constraints and large angle 
radiation. The infra-red cut-off has been used as 



ll\ 
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For ultraviolet cut-off on intermediate gluons, we have 
used [Uj], 
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where /3o = (1 — M 2 jE 2 ) 1 ! 2 and T is temperature of ther- 
mal background. The ultraviolet cut-off on energy for the 
emitted soft gluon has been taken as average momentum 
of the intermediate gluon line as 55( , 



w: 



{91) 



(12) 



Now, the relation between uj and fcj_, uj — kj_ cosh 77, can 
be used to obtain bound on rj from top, which eventu- 
ally excludes all collinear singularities for massless case. 
Finite cut on uj and k± then leads to an inequality, 

coahr)> u> max / kx\ min , (13) 

from which one can easily obtain the bound on 77 as 
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We are now in position to discuss the LPM effect which 
is usually included through a step function 9(ji — r/) 
while evaluating the spectrum of the radiated gluon. It 
basically implies that the formation time of the gluon, 
Tf = {uj) / (k 2 ) must be smaller than the interaction time 
n ~ Aq^ d = 0.49/Tc- This on the other hand imposes 
a restriction on the phase space of the emitted gluon as 
(w) > 2Aq C d ~ 4T C — gT — /j D , provided a s — 0.3, 
T c ~ 170 MeV and the temperature of the plasma, 
T ~ 350 MeV. Thus, the hierarchy in Eq.@ excludes 
the modification of the radiative energy loss due to the 
LPM interference correction through the infrared regula- 
tor, /j£>. Therefore, the present formalism becomes akin 
to the Bethe-Heitler approximation, in which the scat- 
tering centers are well separated and the intensity of the 
induced radiation from different scatterings is additive. 

Now, it is very straightforward to obtain the radiative 
energy-loss through the inelastic processes, viz., Qq(q) — > 
Qq(q)g and Qg — > Qgg [53], for a heavy quark from (|5J|, 
which reads as 
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Equation (|T5j) together with <jTH]) represents radiative en- 
ergy loss of an energetic quark in a canonical way within 
the framework of perturbative QCD along with kinemat- 
ical restrictions for an energetic parton and medium in- 
teraction. 



HEAVY QUARK PRODUCTION IN PP 
COLLISIONS 



At leading order pQCD, heavy quarks in pp collisions 
are mainly produced by fusion of gluons (gg — > QQ) or 
light quarks (qq — > QQ) [56|. The cross-section for the 
production of heavy quarks from pp collisions at leading 
order can be expressed as [56J, y57[ : 



is taken as the baseline for the calculation of the nu- 
clear suppression factor, R AA 44J- Effect of pref actor 



K is diluted during computation of nuclear modifica- 
tion factor due to its identical effects on both initial and 
final distributions profiles. Furthermore, the K-factor, 
if equal for c and b quarks, has not only a diluted ef- 
fect but can actually be neglected in the ratios. The 
shadowing effect is considered using EKS98 parameter- 
ization J58| for nucleon structure functions and here we 
use the CTEQ4M [59j set for nucleon structure function. 
We use Peterson fragmentation function with parameter 
e c = 0.06 and e^ = 0.006 for fragmentation of c quarks 
into D mesons and b quarks into B mesons, respectively. 
All the calculations are done assuming the mean in- 
trinsic transverse momentum of the partons to be zero. 



INITIAL CONDITIONS AND EVOLUTION OF 
THE MEDIUM 

As the heavy quarks are expected to lose most of their 
energy during the earliest time after the formation of 
QGP, we can safely neglect the transverse expansion of 
the plasma while discussing the heavy quark energy loss. 

We consider a heavy quark, which is being produced 
at a point (r, $) in a central collision and moves at an 
angle <f> with respect to r in the transverse plane. If R be 
the radius of the colliding nuclei, the path length covered 
by the heavy quark would vary from to 2R, before it 
exits the QGP. The distance covered by the heavy quark 
inside the plasma in a central collision, L, is given by 
611: 
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where i and j are the interacting partons, f- and / 
are the partonic structure functions and X\ and x 2 are 
the fractional momenta of the interacting hadrons carried 
by the partons i and j. The short range subprocesses for 
the heavy quark production, a — da/dt are defined as: 
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where \A4\ 2 for the processes gg—^QQ and qq^ QQ can 
be obtained from Ref. [561 ] . The running coupling con- 
stant a s at leading order is 
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where Nf — 3 is the number of active flavours and A = 
Aqcd • The pt distribution of production of heavy quarks 
at leading order supplemented with a K-factor sa 2.5 



L(<p, r) = JR. 2 - r 2 sin 
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We can estimate the average distance traveled by the 
heavy quarks in the plasma as: 

Jrdr fL{<t>,r)T AA {r,b = 0)d<t> 

W = ° n° 2 , • ( 21 ) 

Jrdr JT AA {r,b = 0)d<j) 
o o 

where T AA (r, b = 0) is the nuclear overlap function. We 
estimate (L) as 5.78 fm for central Au+Au collisions and 
6.14 fm for central Pb+Pb collisions. 

The temperature of the plasma at a time r, assu ming 
a chemically equilibrated plasma can be expressed as [43 1 



T(r) = 



P(r) 



1.202 (9 N f + 16) 
where the gluon density at time r is given by [43j : 
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Here we consider only the gluon density as the heavy 
quarks lose most of their energy in interaction with glu- 
ons. We also add that the gluon multiplicity is taken as 
3/2 times the number of charged hadrons and the initial 
temperature is obtained using (1221) . assuming an initial 
time. 

We take (^)«1125 for Au+Au collisions at 200 
AGeV H, w 2855 for Pb+Pb collisions at 2.76 ATeV0 
and «4050 for Pb+Pb collisions at 5.5 ATeV \&^. We 
assume that the heavy quark having rapidity in the cen- 
tral region moves along the fluid of identical rapidity. 
This kind o f ap proximation has been used earlier in lit- 
erature [6J-166J. 

We calculate the initial temperature of QGP formation 
T at 200 ATeV as 400 MeV, at 2.76 ATeV as 525 MeV 
and at 5.5 ATeV as 590 MeV, assuming the initial time 
of QGP formation as to = 0.2 fm/c. The critical tem- 
perature T c for the existence of QGP is taken as sa 170 
MeV. The time, by which the plasma will reach the crit- 
ical temperature, r c is found to be sw 2.627 fm/c at 200 
AGeV, 5.9038 fm/c at 2.76 ATeV and 8.375 fm/c at 5.5 
ATeV, assuming Bjorken's cooling law, T 3 r = constant. 

The average path length of the heavy quark inside the 
plasma is calculated as follows. The velocity vt of a 
heavy quark can be expressed as J»j_/ to Tj where my is 
the transverse mass. Thus, the heavy quark would cross 
the plasma in a time Tl = (L)/vt- Now, if t c > tl, 
the heavy quark would remain inside the QGP during the 
entire period, To to tl- But if r c < tl, it would remain 
inside QGP only while covering the distance vt x t c . 
Thus, we further approximate the expanding and cooling 
plasma with one at a temperature of T at r = (L) e ff/2, 
where (L) e s =min [(L), vt X r c ] (see Ref. [43j). 



RESULTS AND DISCUSSION 

In Fig. [2] a comparison of average radiative energy loss 
of an energetic quark traversing in a deconfined quark 
matter produced in Pb-Pb collision at 2. 76 A TeV in the 
present calculation with Djordjevic, Gyulassy, Levai and 
Vitev (DGLV) formalism in Refs. 0, |Hj]. As can be 



seen both light and heavy quarks in the present formal- 
ism, within the gluon emission spectrum of 0(a s ) and 
0(l/k±) as given in (TTJ), lose energy in a similar fashion 
for E > 10 GeV since the effect of mass is small compared 
to the energy. However, it is slightly less than that of a 
light quark for E < 10 GeV, due to the dead cone sup- 
pression at small angles. In addition the results from the 
present calculation differ from that of DGLV 22J, |43j one. 
These differences arise mainly because of the proper kine- 
matic cuts for gluon emission as well as the method used 
to obtain energy loss. The various cuts in the present as 
well as in DGLV formalism are in close proximity except 
the gluon emission in DGLV is constrained only to the 
forward emission angles [27j, 9 < ir/2, whereas in the 
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FIG. 2. (color online): Comparison of average energy loss 
for light quark and charm quark with mass 1.5 GeV in a 
deconfined quark matter produced in Pb-Pb collision at 2.76 
ATeV in the present and DGLV [43] formalisms. For both 
cases the characteristics of the deconfined matter are treated 
in the same footing, i.e., the strong coupling a s — 0.3 and the 
average path length, (L) ~ 6.14 fm, traversed by an energetic 
quark in a deconfined medium produced in such collisions. 
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FIG. 3. (color online): Same as Fig. [2] but only for charm 
quark in Au-Au collision at 200 AGeV with (L) - 5.78 fm. 



present calculation the full range of 6 is taken care off 
through the variable ij as shown in (fT3|) and (|14[) . 

In Figs 01 and 2] we have displayed average energy loss 
of a charm quark in a deconfined quark matter, respec- 
tively, at 200 AGeV Au-Au collision at RHIC and 5.5 
ATeV Pb-Pb collision at LHC. We find that at RHIC 
energies the average energy loss of a charm quark in our 
formalism is higher than that of the DGLV formalism for 
the considered energy range, (0 < E < 50) GeV, of the 
charm quark. On the other hand FigU is qualitatively 
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FIG. 4. (color online): Same as Fig. [3] in Pb-Pb collision at 
5.5 ATeV with (L) = 6.14 fm. 



FIG. 6. (color online): Nuclear modification factor Raa for 
D mesons with both collisional and radiative energy loss in 
Pb-Pb collision at 2.76 ATeV. The data are from Ref. H3] but 
only the systematic error bars are shown here. 
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FIG. 5. (color online): Collisional energy loss of charm 
quark 31j in Pb-Pb collision at 2.76 ATeV and 5.5 ATeV 
at LHC, and 200 AGeV at RHIC energies. 
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FIG. 7. (color online): Nuclear modification factor, Raa, for 
D mesons in Pb-Pb collision at 5.5 ATeV. 



similar to Fig{2] in terms of comparison of two formal- 
ism for heavy quark. As seen the average energy loss of 
charm quark is larger in the present formalism only in 
the domain, (0 < E < 15) GeV, of the charm quark and 
beyond which it is less compared to the DGLV formal- 
ism. The difference, in fact, increases as energy of the 
quark increases. 

In Fig. [5] we display a comparison of collisional energy 
loss of charm quark as calculated by Peigne and Peshier 
(PP) in Ref. [3l| for RHIC and LHC energies. As seen 
the collisional energy loss increases with the increase in 
centre of mass energy of the colliding ions. 

In Fig. [5] the nuclear suppression factor, Raa, for 
D meson is displayed considering both radiative and 



collisional energy loss and compared with the ALICE 
data 67( at 2.76 ATeV. As can be seen the differences 
in radiative energy loss between the present and DGLV 
formalism discussed in Fig. [2] for 2.76 ATeV in Pb-Pb 
collisions is clearly reflected in Fig. [6] For the present 
calculation it is manifested in gradual increase of Raa 
of D meson 67j for transverse momentum, p± > 5 GeV 
whereas in DGLV case it remains almost constant. The 
suppression factor obtained in the present formalism with 
radiative energy loss is in close agreement with the most 
recent data from ALICE collaboration at 2.76 ATeV J67J . 
On the other hand the inclusion of the collision contri- 
bution is found to suppress Raa further in both cases. 
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FIG. 8. (color online): Raa with only radiative energy loss 
for non-photonic single electron from the decay of individual 
D mesons and B mesons in Au-Au collision at 200 AGeV. 
The data are from Ref. [g|. Both systematic and statistical 
error bars are shown for STAR data whereas only systematic 
error bars are displayed for PHENIX data. 



FIG. 9. (color online): Raa with collisional and radiative 
energy-loss for non-photonic single electron from the com- 
bined decay of both D and B mesons in Au-Au collision at 
200 AGeV. The data are from Ref. Q. Both systematic and 
statistical error bars are shown for STAR data whereas only 
systematic error bars are displayed for PHENIX data. 



As found the data suggest that the collisional contribu- 
tion may be small. Nonetheless, more data in the high 
p± domain is necessary to know the actual trend of the 
energy loss of charm quark and will finally constrain the 
various energy loss and jet quenching model in the liter- 
ature. We also expect a similar rise in light hadrons for 
high p± since both light and heavy quark lose energy in a 
similar fashion as shown in Fig.[2j However, we note that 
the ALICE data on Raa for inclusive charge hadrons [68( 
at 2.76 ATeV in Pb-Pb collision has also shown a similar 
increasing trend as p± increases. It is natural to believe 
that such data is completely dominated by the contri- 
bution from light hadrons. For completeness, we also 
display Raa for LHC energy at 5.5 ATeV in Fig. [7] 

In Fig. |8]the nuclear suppression factors for individual 
decay of D and B mesons to non-photonic single electron 
is displayed considering only the radiative energy loss for 
RHIC energy at 200 AGeV. As expected the contribution 
from the B decay is small compared to that of D decay. 
In Fig. |9] the total contribution of single electron from 
D and B decay is shown considering both radiative and 
collisional energy loss. It is found that the contributions 
of the collisional energy loss is important at RHIC en- 
ergy. We also compare our results with that of DGLV. In 
Fig. 1 101 we give prediction for single electron result for 
LHC energy at 2.76 ATeV. 



SUMMARY 

We obtain the radiative energy loss of a heavy quark 
akin to the Bethe-Heitler approximation by considering 
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FIG. 10. (color online): Same as Fig. |9]in Pb-Pb collision at 
2.76 ATeV. 



the most generalised gluon emission multiplicity expres- 
sion derived very recently. This suggests that both en- 
ergetic heavy and light quark lose energy due to gluon 
emission almost similarly and the mass plays a role only 
when the energy of the quark is of the order of it. The 
hierarchy used for simplifying the matrix element as well 
as for obtaining the gluon radiation spectrum imposes 
a restriction on the phase space of the emitted gluon in 
which the formation time is estimated to be less than 
the interaction time. This suggests that the LPM inter- 
ference correction may be marginal. Further, we compare 
our results with the DGLV formalism and it is found to 



differ significantiy. To compute the nuciear suppression 
factor for Z?-meson we consider both radiative and colli- 
sion energy loss along with longitudinal expansion of the 
medium. The nuclear modification factor for D-meson 
with radiative energy loss obtained in the present formal- 
ism has an increasing trend at high p± and found to agree 
closely with the very recent data from ALICE collabora- 
tion at 2.76 ATeV. When the collisional counter part is 
added independently, the further suppression is obtained 
in the nuclear modification factor. This suggest The non- 
photonic single electron data at 200 AGeV RHIC energy 
requires contributions from collisional energy loss as well 
from B decay. However, it is necessary to obtain both 
radiative and collisional energy loss from the same for- 
malism to minimize the various uncertainties, which is 
indeed a difficult task. Moreover, data at high p± region 
with improved statistics are required to remove prejudice 
on different energy loss and jet quenching models. 
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